Classical Spin Models with Broken Continuous Symmetry: Random Field Induced 
Order and Persistence of Spontaneous Magnetization 
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We consider a classical spin model, of two-dimensional spins, with continuous symmetry, and 
investigate the effect of a symmetry breaking unidirectional quenched disorder on the magnetization 
of the system. We work in the mean field regime. We show, by numerical simulations and by 
perturbative calculations in the low as well as in the high temperature limits, that although the 
continuous symmetry of the magnetization is lost, the system still magnetizes, albeit with a lower 
value as compared to the case without disorder. The critical temperature at which the system 
starts magnetizing, also decreases with the introduction of disorder. However, with the introduction 
of an additional constant magnetic field, the component of magnetization in the direction that is 
transverse to the disorder field increases with the introduction of the quenched disorder. We discuss 
the same effects also for three-dimensional spins. 



I. INTRODUCTION 

Disordered systems, both classical and quantum, lie at 
the center-stage of condensed matter physicaii^. Chal- 
lenging open questions in disordered systems include 
those in the realms of spin glasses^, neural networks^, 
percolation^, and high superconductivity^. Phenom- 
ena like Anderson localization^!^, and absence of magne- 
tization in several classical spin models^'^S'ii are effects 
of disorder. 

In particular, classical ferromagnetic spin models with 
discrete, or continuous, symmetries are very sensitive to 
random magnetic fields, distributed in accordance with 
the symmetry, in low dimensions^. For instance, an arbi- 
trary small random magnetic field with Z2 (±) symmetry 
destroys spontaneous magnetization in the Ising model in 
2D at any temperature T, including T — Q. Similar effect 
holds for the XY model in 2D at T = in a random field 
with t/(l) {SO{2)) symmetry, or Heisenberg model in 2D 
at T = in SU (2) (S'0(3)) symmetry in random field. In 
these cases, the effects of disorder amplify the effects of 
continuous symmetry, that destroys spontaneous magne- 
tization at any T > 0. The effect is even more dramatic in 
3D, where the random field destroys spontaneous mag- 
netization at any T > 0. (For a general description of 
these, seeSiiaiii.) 

The appropriate symmetry of the random field is essen- 
tial for the above mentioned results. The natural ques- 
tion arises as to what happens if the distribution of the 
random field does not exhibit the symmetry, in particular 
the continuous symmetry. Yet another natural question 
is how does the spin systems in random fields behave in 
the quantum limit. The latter question is particularly in- 
teresting in view of the fact that nowadays it is possible 
to realize practically ideal models of quantum spin sys- 
tems (with spin s=l/2, 1, 3/2, and with Ising, XY, or 
Heisenberg interactions) in controlled random fields^ii^. 
It is therefore very important to understand the physics 
of both classical and quantum spin models in random 



fields that break their symmetry. 

In this paper, we will consider the classical XY spin 
model in a random field that breaks the continuous U(l) 
(50(2)) symmetry. We investigate this model in the 
mean field approximation^'^. Despite its simplicity, this 
model with the two-dimensional spin variable, magne- 
tizes in the absence of disorder below a certain criti- 
cal temperature, which can be calculated exactly. As 
a result of continuous symmetry, the possible values of 
spontaneous magnetization form a circle in the plane. A 
random magnetic field pointing along the y-direction is 
introduced, by adding a new term to the energy of the 
model. This term breaks the continuous symmetry of 
the model, but the critical temperature persists. The 
present paper studies the critical behaviour and proper- 
ties of spontaneous magnetization in the resulting mean- 
field disordered system. We prove that, as may be ex- 
pected, adding a random field lowers the critical tem- 
perature. Next, we show that the magnetization of the 
disordered system is lower than that of the pure one — 
another intuitively plausible result, first shown numeri- 
cally and then by perturbation expansions: one around 
the critical temperature of the pure model and another 
for low temperatures. 

Next, we introduce a constant magnetic field, which 
breaks the continuous symmetry of the model even in 
the absence of disorder. In fact, the system now mag- 
netizes at all temperatures in the direction parallel to 
the magnetic field. When we now also add a random 
field as described above, the length of the magnetiza- 
tion vector decreases again. Moreover, the magnetiza- 
tion gets atrracted towards the X- axis, i.e. the direction 
transverse to that of the random field. However, the X- 
component of the magnetization can increase for certain 
choices of the constant field. We view this effect as a case 
of "random field induced order" , by analogy with the ef- 
fect studied in^^, where numerical evidence was given for 
appearance of magnetization in the XY model on a two- 
dimensional lattice with the introduction of disorder. In 
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contrast to the present work, in this other case no mean- 
field approximation was used and no uniform magnetic 
field was introduced. 

The effect of "random field induced order" has, of 
course, a long historj*i^. Recently it has become vividly 
discussed in the context of XY ordering in a graphene 
quantum Hall ferromagnet"*^^, and ordering in '^He-A 
arerogel and amorphous ferromagnets^^. Let us stress 
that the novel aspects of our previous paper— consists 
in clarifying certain aspect of the rigorous proof of the 
appearance of magnetization in XY model at T = 0, 
presentation of a novel evidence for the same effect at 
T > 0, and a proposal for realisation of quantum ver- 
sion of the effect with ultracold atoms. In the subse- 
quent paper—, we have shown how the "random field in- 
duced order" exhibits itself in a system of two-component 
trapped Bose-Einstein condensate with random Raman 
inter-component coupling. The novelty of the present 
paper lies in systematic mean field treatment of the dis- 
ordered XY model with particular emphasis on the re- 
sponse to the constant magnetic field. 

The paper is arranged as follows. In Sect. lUthe ferro- 
magnetic XY model is introduced. A symmetry breaking 
random field is added in Sect. IIIII and the results of nu- 
merical simulations of and perturbative calculations on 
the resulting model are presented. In Sect. llVi study the 
system with an additional constant field and, in partic- 
ular, show presence of random field-induced order. We 
discuss the results in Sect. |Vl arguing in particular that 
the analogs of our results will hold for the mean-field 
version of the classical Heisenberg model. 



II. FERROMAGNETIC XY MODEL: 
FIELD APPROACH 



MEAN 



Consider a lattice, each site i of which is occupied 
by a "spin", which is a unit vector cfi = (cos 0^, sin 0^) 
on a two-dimensional plane (called the XY plane). The 
nearest-neighbor ferromagnetic XY model is defined by 
the Hamiltonian 



H 



XY 



E 



(1) 



with a coupling constant J > 0. This model 
does not have any spontaneous magnetization, at any 
temperature, in one and two dimensions (Mermin- 
Wagner-Hohenberg theorenii^), while a finite magneti- 
zation appears in higher dimensions for sufficiently low 
temperaturesiSiSS. 

Let us assume that the total number of spins in our 
system is N. In the mean field approximation every spin 
is assumed to interact with all other spins (not just with 
the nearest neighbors) with the same coupling constant 
— J. Therefore, the contribution of the spin at i to the 



total energy of the system equals 

where we divided the energy term by iV in order to pre- 
serve its order of magnitude. This effective interaction, 
replacing the nearest neighbor interaction in Hxy , is for 
large N approximately equal 



(2) 



where fh — '^i- ^^'^ mean field approximation 

consists of treating m as a genuine constant vector and 
adjusting it so, that the canonical average of the spin at 
(any) site i equals this constant. If the system is in canon- 
ical equilibrium at temperature T, the average value of 
the spin vector Oi is 



J ex:p{(3Jrh ■ (Ji)d<Ti 



(3) 



where /3 = l/^ksT), with ks being the Boltzmann con- 
stant. This average is independent of the site i. Con- 
sistency requires that the left hand side (l.h.s.) of the 
above equation be equal to the magnetization m. Hence, 
we obtain the mean field equation 



/ tr exp(/3 Jm • (7)d(T 
J exp(/3Jm • (7)dd ' 



(4) 



where we have dropped the index i. Equations of this 
type, for various modifications of the original interaction 
are the main subject of this work. 
Let 



m — (m cos a, m sin a) . 



(5) 



For sufficiently high temperatures, the only solution fo 
the mean field equation is to = 0. There exists a /3", 
such that for /3 > this system magnetizes. (Later, 
we will consider the case of a system with an additional 
quenched disordered field of strength e. The superscript 
of is anticipation of that case.) By symmetry, the 
solutions of the above mean field equation (Eq. ^) form 
a circle 



mo, 



(6) 



with a strictly positive radius toq, for any (3 > (3^. Choos- 
ing the phase a = 0, Eq. ([4]) reduces to 



cos 9 exp{pjm cos 
exp{(3Jm cos 9)d9 



(7) 
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FIG. 1: (Color online.) The behavior of magnetization (ver- 
tical axis) with respect to /3J (horizontal axis). The (red) 
squares represent the magnetization of the XY model without 
disorder (the corresponding Hamiltonian being Hxy), while 
the (blue) traingles are that for the same, but with disorder 
(the corresponding Hamiltonian being i/e with e = 0.1). 



A. Quenched averaging 

The disorder considered in this paper is "quenched", 
i.e. its configuration remains unchanged for a time that 
is much larger than the duration of the dynamics consid- 
ered. In the systems that we study, it is the local mag- 
netic fields that are disordered. They are random vari- 
ables following certain probability distributions. Since 
the disorder is quenched, a particular realization of all 
the random variables remains fixed for the whole time 
necessary for the system to equilibrate. An average of a 
physical quantity, say is thus to be carried out in the 
following order. 

(a) Compute the value of the physical quantity A, with 
the fixed configuration of the disorder. 

(b) Average over the disordered parameters. 

This mode of averaging is called "quenched" averaging. 
It may be mentioned that an averaging in which items 
(a) and (b) are interchanged in order, is called "annealed" 
averaging. Physically it corresponds to a situation when 
the disorder fluctuates on time scales comparable to the 
system's thermal fluctuations. 



where we have taken 

(7= (cos 6*, sin 0). (8) 

In Fig. [2 the red squares represent the cross- 
section, of the surface of solutions of Eq. ([4]) in the 
(to cos a, m sin a, /3) Cartesian space, in the cos a ~ 
plane. 

From numerical simulations (see Fig. [T]), we found 
that J K, 2.00. One can show analytically that J is 
exactly 2, as follows. Let us denote the right hand side 
(r.h.s.) ofEq. F{m). The condition for the system 

to magnetize, for a given value of /?, is that the derivative 
of the r.h.s. of Eq. ([7]) at to = 0, i.e. F {m)\m=o should 
be greater than the derivative of the l.h.s. It is easy to 
check that 

F'{m)\ra=a^PJ/2. 

This implies that the system possesses a nonzero magne- 
tization, if and only when 

/3 > 2/ J = Pi 

III. FERROMAGNETIC XY MODEL IN A 
RANDOM FIELD 

We will now consider the effect of additional quenched 
random fields. Let us begin by reminding the notions of 
quenched disorder and quenched averaging. 



B. The model and the mean field equation for 
magnetization 

The XY model with an inhomogeneous magnetic field 
has the interaction: 

H = - J ^ (Ti ■ (fj - e^^hi ■ di. (9) 

where the two-dimensional vectors hi are the external 
magnetic fields, up to a coefficient e. In the sequel hi are 
random variables of order one, they model the disorder in 
the system and thus e measures the disorder's strength. 
More precisely, let hi be independent and identically dis- 
tributed random variables (vector- valued) . We want to 
study the effect of including such a random field term in 
the XY hamiltonian at small values of e. As argued irti^, 
in lattice XY models this effect depends critically on the 
properties of the probability distribution of the random 
fields. 

If the distribution of the hi is invariant under rotations, 
there is no spontaneous magnetization at any nonzero 
temperature in any dimension d < 4^ ii0iii 

We now want to see the effect of a random field that 
does not have the rotational symmetry of the XY model 
interaction ([1]), considering the case when 

hi = ni ■ Sy (10) 

where rji are scalar random variables with a distribu- 
tion symmetric about and ey denotes the unit vector 
in the y direction. The main result ofl^ is that on the 
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two-dimensional lattice such a random field will break 
the continuous symmetry and the system will magnetize, 
even in two dimensions, thus destroying the Mermin- 
Wagncr-Hohcnberg effect. Above two dimensions the 
pure XY model magnetizes at low temperatures and it 
has been suggested in^ that the uniaxial random field 
as described above may enhance this magnetization. In 
the present paper we want to study related effects at the 
level of a simpler, mean-field model, which allows for a 
more detailed analysis and more accurate simulations. 
We consider the mean-field Hamiltonian given by 



where, as ffi before, 



-Jm ■ a — eff ■ a, 



(11) 



is the quenched random field in the y-dircction. Here rj 
is a scalar, symmetric random variable, which we assume 
here to be Gaussian distributed with zero mean and unit 
variance, e (> 0) is the parameter (typically small) that 
quantifies the strength of the randomness. 

The corresponding mean field equation for magnetiza- 
tion is: 



= AVr, 



J (Texp(/3JTO • d + Pe7]ay)dd 



J exp{(3Jrh ■ a + (3triay)d(T 



(12) 



Here Avr^{-) denotes the average over the disorder, i.e. 
the integral over rj with the appropriate distribution (here 
assumed to be unit normal). 



C. Numerical simulations 

It follows from the symmetry of the distribution of t] 
that all solutions of the equation (Eq. (H^)) have zero 
y-component. In the case of e ^ 0, the system again 
does not magnetize at high temperature, as in the case 
of e = 0. However, there exists a critical temperature, be- 
low which a transverse (with respect to the direction of 
the random field) magnetization appears. More precisely, 
there exists a fil such that for P > (3^, the magnetiza- 
tion equation has two solutions with zero y-components, 
whose X-components equal m and — m (where m > 0) 
(along with the trivial solution m — 0). We will study the 
dependence of m on the temperature and on the disor- 
der strength e. Therefore, in contrast to the continuous 
set (a circle) of solutions of the mean field equation in 
the system without disorder (modelled by Hxy), in the 
disordered case there are just two possible values of the 
spontaneous magnetization. In Fig. [1] the blue trian- 
gles represent the magnetization for e = 0.1, while the 
red squares correspond to e = 0. The same color code 
applies to Fig. [2l where, in addition, pink stars repre- 
sent the case e = 0.15, and green circles — e = 0.2. All 
the curves show two real solutions (m and — m), of the 
corresponding mean field equation, at low temperatures. 
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FIG. 2: (Color online.) The magnetization (m) with respect 
to /3J (horizontal axes): the (red) squares, (blue) triangles, 
(pink) stars, and (green) circles are for the system described 
by the Hamiltonian with e = 0, e = 0.1, e = 0.15, and 
e = 0.2 respectively. The figure on the left correspond to the 
behavior of the spontaneous magnetization near the critical 
temperature, while that on the right correspond to the same 
at lower temperatures. The decrease of magnitude of the 
magnetization due to disorder, is clearly of the order of e^. 



From Figs. [l]and[2l it is clear that m{(3, e) is smaller than 
Too (see Eq. (j6|)), at low temperatures. They coincide for 
high temperatures, as both of them are vanishing in that 
regime. Numerical simulations show that the difference 
S{f}, e) — Too — to(/3, e) — to(/3, 0) — m{(3, e) > is of the 
order of e^, in the regime of /3 > (see Fig. 



D. Scaling of critical temperature and 
magnetization with disorder: Perturbative approach 

In this subsection, we will study the mean field Hamil- 
tonian in Eq. (jlip . using perturbation theory and com- 
pare the results with the numerical results of the last 
section. Since, as argued earlier, the spontaneous mag- 
netization can only have nonzero X-component at any 
temperature, the mean field equation (Eq. (I12p ) reduces 
to 



AVr 



cos 9 exp{pjm cos 9 + Per] sin 9)d9 
J^^ exp{l3Jm cos 9 + Per] sin 



= F,(to).(13) 



Note that F^{m) = F(to) for e = 0. 
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1. Critical temperature 

To find the critical temperature, a similar method as 
in the case of e = (in the preceding section) is applied. 
The condition for non-zero magnetization is now given 

by 

P'M-^-^+o{e')>l, (14) 

where o(C") denotes a term which is of order higher than 
(^". This impHes 



(15) 



Therefore, we obtain negative corrections to the criti- 
cal temperature, as observed in the numerical simulations 
(see Fig. [2]). 



2. Scaling of magnetization near criticality 

The magnetization goes to zero as the temperature ap- 
proaches the critical temperature. Let us now use pertur- 
bation techniques to see the behavior of magnetization m 
near criticality. 

Before considering the disordered case (e ^ 0), let us 
first consider the case when e = 0. In this case, we expand 
F{m) (r.h.s. of Eq. ([7])) in m around to = 0: 

'I 

I Til III 

F(m) = mF (0)-| F (0)-f higher order terms. (16) 

6 

Differentiating Eq. ((4]) and considering the resulting ele- 
mentary integrals one can easily see that, by symmetry, 
F(0) = F"{0) = and that 



F (0) = /3J/2, F (0) = -3/3V78 



(17) 



Putting these values in Eq. (|16p. we obtain that the 
magnetization near f3^ equals 



Too 



2V2 



(18) 



plus higher order terms. 

A similar technique can now be used to study the be- 
havior of TO near f3c in the presence of disorder. Again 
F^{Q) = F^ (0) = 0. Hence the expansion of F^{m) near 
zero will be 



F,(m) = toFJO)- 
where 

F'M 

f:"{o) 



-F^ (0)+higher order terms, (19) 
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Putting these derivatives in Eq. (fT9|) . we obtain the cor- 
rection of magnetization due to disorder as 



This is in full agreement with numerical simulations, that 
also showed a decrease of magnetization of order e^, in 
the disordered case, as compared to the case when e = 0. 



3. A modified Bessel function and its expansion for large 
arguments 

In the following, we will have numerous occasions to 
use the modified Bessel function 



Iniz) 



exp(z cos d) cos{nd)d^, 



(20) 



where n is an integer, and z is a real number, and the 
expansion of /„(z) for large \z\i^ 



In{z) 



exp(z) r ^ - 1 (pi - 1)(m- 9) 



'2tiz 



8z 



2!(8z) 



(M~l)(M-9)(/i-25) 
3!(8z)3 



o({l/zy 



(21) 



where n is fixed and /i = 4n^. Actually, the function and 
its expansion are true^ for certain complex ranges of the 
parameter z. However, we will only use them for real z. 
Here o{{l/v) ) denotes an expression, containing terms 
of order higher than (i)": (1/f )"+^ , (l/i;)"+^ , . . .. 



4- Scaling of magnetization at low temperature 

We now study behavior of to at low temperatures, i.e. 
for large /?. 

We again start from the case e = 0. Note that the 
numerator and denominator of F{m) (after some simple 
modifications) are of the form of In{z). Therefore, for 
large /3, we can use the asymptotics of the Bessel func- 
tion in Eq. Ij2ip to obtain a low temperature expansion 
of F{m). Using this expansion, we obtain the following 
equation for to, from Eq. ([7]): 



m^-m^ + — +o{l/P) = Q. 



Since m 1 as /3 — > cxd, let us write m as 
m=l-^+o[l/P). 



(22) 



(23) 



Putting this in Eq. ([22|) . we finally obtain the behavior 
of the magnetization for the case when e = 0, for large (3: 



Too = 1 



1 

2J^ 



o(l//3). 



(24) 
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Using the same technique for the disordered case, we 
expand the numerator and denominator of F^{m), using 
Eq. ([21]), obtaining 



1 



e 

2J2 



1 

2J 



+ t 



1 

2J2 



mo 



1^- 


1 


\(—-—)] 


l2J2 







ho(l//3) 

-o(l//3). 
(25) 



As before, disorder leads to corrections of order to 
the magnetization at low temperature. This is again in 
agreement with the numerical simulations (see Fig. [2]). 



IV. FERROMAGNETIC XY MODEL IN A 
RANDOM FIELD PLUS A CONSTANT FIELD: 
RANDOM FIELD INDUCED ORDER 

We have seen in the preceding section that a random 
field that breaks the symmetry of the XY model, restricts 
possible magnetization values to a discrete set. Although 
the system still magnetizes, we no longer have continu- 
ously symmetry of the set of solutions to the mean field 
equation, as a result of adding a symmetry-breaking ran- 
dom field. In this section we explore the effects of such 
a random field on a system which already has a unique 
direction of the magnetization, determined by a uniform 
magnetic field. 

First consider the case in which the planar symmetry in 
the XY model is broken by applying a constant magnetic 
field h alone. That is, according to the general mean field 
strategy, we are looking for the solutions of the following 
equation: 



/ (Texp(/?JTO ■ a + l3h ■ (T)da 
J e'}q){f3Jrn ■ a + (3h ■ a)du 



(26) 



Let h = (ft, cos X, hsvax). We suppose that < ft. < 1 and 
Q < X < 7r/2. As expected, due to the applied constant 
field, the mean field equation has a unique solution of the 
magnetization m at all temperatures, and the solution 
is a (positive) multiple of ft, but of reduced magnitude. 
Red squares in Figs. O [21 El and [H correspond to the 
magnitude (m) and the cosine of the phase (a) of the 
magnetization m. 

Let us now, in addition, apply a random field {77} in 
the Y-direction. The new mean field equation is 



m 



AVn 



J (jexp(/3Jm ■ d + ph ■ a + (3triay)da 
J e:>q){l3Jm ■ a + l3h ■ a + I3erjay)da 



(27) 

Here we have to solve the two simultaneous equations, 
given by Eq. ([27)) . to obtain the magnitude and the 
phase of the magnetization vector m. In all the previous 
mean field (vector) equations, one could apriori predict 
the phase of the magnetization. Just as in the case of 
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FIG. 3: (Color online.) Plot of cos(a) with respect to /3J. Red 
squares represent the case when the XY model has the applied 
constant field h with h = J and x = n/4. Blue triangles 
are for the same system but with the additional symmetry 
breaking random field of strength e = O.IJ. 
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FIG. 4: (Color online.) This is the same plot as in Fig. [3] 
but for X = tt/6. 



a constant field h and e = 0, again the solution remains 
unique. 

Just as in the previous sections, we will now compare 
the magnetization of the system without disorder (i.e. 
e = 0, and for which the mean field equation is given 
by Eq. ([^ )). with the system in which e ^ (and for 
which the mean field equation is given by Eq. ([27[) ) keeping 
h strictly positive in both cases. Let us denote the two 
Hamiltonians by Hh and Hh.e respectively. We do the 
comparison by numerical simulations as well as pertur- 
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FIG. 7: (Color online.) Schematic diagram of the magneti- 
zation of XY ferromagnets without and with disorder, in the 
presence of a constant magnetic field. The figure on the left 
indicates the behavior of m in the presence of h, but when 
e = 0, while the one on the right is when there is a positive e. 



FIG. 5: (Color online.) Plot of the magnitude m of the mag- 
netization with respect to /3J. Red squares represent the case 
when the XY model has the applied constant field h with 
h = J and x — 7r/4. Blue triangles are for the same system 
but with the additional symmetry breaking random field of 
strength e — O.IJ. 
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FIG. 6; (Color online.) This is the same plot as in Fig. [5] 
but for X = 7r/6. 



scribed by Hh). This is seen from numerical simulations 
(see Figs. [5]and[6]), as well as by perturbation techniques 
at low temperatures. In addition, numerical simulations 
(as shown in Figs. [3] and 2]) show that the cosine of 
the phase of the magnetization, i.e. cos(a) increases in 
presence of the random field. Therefore, the phase of 
the magnetization vector moves towards the X-direction 
(i.e. the direction transverse to the applied random field). 
This is also corroborated by perturbative approach at low 
temperatures. The schematic diagram in Fig. [7] shows 
the change of behavior of the length and phase of the 
magnetization with and without disorder, in the presence 
of a constant field. 

The y-component, my — to sin a, of the magnetization 
has the same relative behavior as the length m, in sys- 
tems described by Hh and i?/i,£, i.e. for small e > it is 
lower than for e = 0. 

However, the X-component, rrix — to cos a, of the mag- 
netization, 771, behaves in a very interesting way. Its value 
in the system described by Hh^e can be both higher as 
well as lower than its value in the system described by 
Hh- The numerical simulations in this case are given in 
Fig. E 



batively at low temperatures (Sect. IIVAI below). A per- 
turbation approach, similar to the one in Sect. IIII D 21 
can be done at high temperatures also. We refrain from 
doing it, as the high temperature behavior in this case is 
not so interesting, in view of absence of a phase transi- 
tion. 

The length m of the magnetization vector is shrunk 
inthe system described by Hfi ,:, compared to the case in 
which there is no disorder in the system (i.e. the one de- 



A. Magnetization at low temperature: 
Perturbative approach 

To obtain the behavior of magnetization at low temper- 
ature, we will use the implicit function theorem, which 
we now state. Let an equation f{xi,X2) = of two 
variables xi and X2 be such that f{xi,X2) = at 
(x 1,2:2) = {x\,X2)- X2 is in general an unknown func- 
tion of xi. But we may still understand the character of 
^^\xi=x1 : by using the fact that (under certain regularity 
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where all the total and partial derivatives are taken at 
e = 0. The above system of equations can be solved for 

the second (total) derivatives and -3^, at e = 0, 

once we can find the partial derivatives at e = 0. 

The partial derivatives in Eq. are calculated using 
the following strategy. We have 



(33) 



J/k^T 



FIG. 8: (Color online.) Order from disorder: Red squares 
are the X-components of the magnetization in the absence 
of disorder, while blue triangles are those in the presence of 
disorder e — O.IJ. The constant field is present in both cases, 
with X = Tv/G. 



conditions on / near {xi,X2)) 



dl 

dxi 



5/ 

0X2 



dx2 

{x1,x°) dxi 



(28) 



The usual statement of the implicit function theorem is 
that when is nonzero at {xi,X2)), we can solve the 
equation f(xi,X2) = for X2 uniquely near this point 
and the derivative of the resulting function (x2 as a func- 
tion of xi) at Xi can then be calculated from the above 
equation. However, in the case when the first derivatives 
vanish at a certain point, we can use a simple extension 
of it to calculate the second derivatives. Such a situation 
appears in the calculations below of the second deriva- 
tives of the magnetization with respect to e. 

The mean field equations that we work with here can 
be written in the form 



where 



with 



^ rri. — 



d 



d 



dm-r dr 



(29) 



(30) 



r = logg y" exp (-/3i/,i) or log^ J exp{-f3Hh,e) ■ 

(31) 

It follows from symmetry of the distribution of rj that 
rh is an even function of e and, consequently and 



vanish at e = 0. 



It follows that 



de2 



1 - 



1 d^T^ 



(3 J dm^ 
1 d^T 



[3 J dm^ 



1 

]3J 
1 

]3J 



dh 



dmydrrix de^ 



_dh_ 



d'^edrriy dniydnix dt^ 



where for any observable A, {A) is the Gibbs average 
/ Aexp(-/3H) 



{A) 



(34) 



with H being the relevant Hamiltonian {H^ or H^ ,:). Of 
course, in the case of the system described by the Hamil- 
tonian Hfi, the quenched averaging with respect to i] is 
not required. Using this notation we have, differentiating 
the formula for T twice. 



1 d dT 

(3 J de dnix 

and 

1 92 QY 



Avr, [p-n ( (cos e sin 9) - (cos 0) (sin 0))] , 

(35) 



AVr 



(3 J de^ dm, 

2 (cos ^ sin 0) (sin 0) -f 2 (cos 0) (sin 9) ^ - (cos 9) (sin^ 0) 



(36) 



We expand these partial derivatives with respect to 

at e = 0, using the expansion of the modified Bessel 

function /„. After some calculations, we obtain 



and 



dm,. 



de2 



= (.,^)+ order of i, (38) 



where the functions X and Y are given by (for j — J /h) 

AD + BC 



X{x,j) 
Y{x,j) 



BD + AC 
DE-C^' 



(39) 
(40) 



where 



(32) 
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1 r 1 5 1 

— ((37 + 1) cos a; + (457 + 63) cos 3a:) (j + 3) sin2xsinx H — (i + 2) cosxcos2x 

(j + 1)3 L 32 8 4 

-i(3j - 7)cosa;sin^x , (41) 
8 J 

1 3 3 1 

- -^(3(3^ + l)sina; + (45j + 63) sin 3a;) + -(j + 2) cos2a;sina; + -(j + 3) sin^ a; , (42) 



(j + 1)^ 




FIG. 9: (Color online.) Plot of the function X{x,j) with 
respect to x and j = J/h. Note that there are ranges of the 
{x,j), for which the function X is positive. This fact gives 
rise to the phenomenon of random field induced order in the 
system described by the Hamiltonian Hh.t- 



FIG. 10; (Color online.) Plot of the function Y{x,j) with 
respect to x and j. Clearly, it is negative for the entire range 
of X and j. 



C{x,j) = - 



J cos X sm X 





- 1 




j cos^ 


a; 4 


^ 1 




- 1 




j sin^ 


a; 4 


- 1 



J + 1 



(43) 
(44) 
(45) 



We plot the functions X and Y in Figs. [9landfT0l 

Therefore, at low temperatures, we have, up to order 

= TOa;|e=o + \T2^ ( ^' T ) + 4 ) ' (46) 



,|.=o + ( U + order of i ) . (47) 



From Fig. [TOl it is clear that the y-component of the 
magnetization always decreases in the presence of disor- 
der. However, Fig. [5] shows that there are ranges in the 
parameter space (a;, j), for which the quenched averaged 
X-component, rux, of the magnetization increases in the 
presence of disorder, compared to the case when there is 
no disorder. As noted before, this is in agreement with 
our numerical simulations. 

We have also considered the effect of disorder on the 
length m and phase a of the magnetization. For the 
phase, we consider the expansion of tan(a) = which 
is given by 



tan(a) — — - 
mx 



+ o(e2), (48) 
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FIG. 11: (Color online.) Plot of the function S{x,j) with 
respect to x and j. It is again negative for the entire range of 
X and j. This is in agreement with the numerical results in 
Figs. OandU 



FIG. 12: (Color online.) Plot of the function [Xm^ + 
Ymy)\^^o with respect to x and j. It is negative for all x 
and j, and numerical results in Figs. [5]and|6]corroborate this 
effect. 



with 



DISCUSSION 



d my d rrij: 



e=0 



1 1 



2 1 .2"5'(a;,j) + order of -, 



where 



S{x, j) — Y{x, j) cos X — X{x, j) sin x, 



(49) 



(50) 



with X and Y being given by Eqs. ((23) and ((M)) . 

As shown in Fig. [TI] S (x, j) is negative for all x and 
j. Consequently, the phase a always bends towards the 
X-direction in the presence of disorder (since tan(a) de- 
creases in a, cos (a) increases), as we have already seen in 
simulations. Note that < a < 7r/2. The square of the 
length of the magnetization is given by (up to order e^) 



2 I 2 / 2 I 2m 

+2e^ ( {Xrux + Ytriy^^Q + order of 



(51) 



As seen on Fig. [TH {Xmx + Ymy)\^=a is always negative, 
showing that the length of the magnetization decreases 
in the presence of disorder. Note that the behavior of the 
length and phase obtained perturbatively, matches what 
is shown schematically in Fig. [T] 



To summarize, we have considered classical systems of 
two dimensional spins, and studied the interplay between 
continuous symmetry and symmetry-breaking quenched 
disordered field, in the mean field approximation. We 
found that in case of a system in a uniform magnetic 
field, disorder may enhance one component of the order 
parameter. 

In this paper, we have explicitly considered only the 
situation when the spins are two-dimensional. It is natu- 
ral to ask analogous questions for three-dimensional spins 
with continuous symmetry. Below we argue that for the 
canonical system of this kind — the classical Heisenberg 
model — the behavior is similar to that of the XY model. 

To study the behavior of magnetization in the lattice 
Heisenberg model in the presence of disorder, we put at 
all sites random fields in the Y-direction. The Hamilto- 
nian of the resulting disordered system is given by 



H 



H 



(52) 



where ff = rj ■ Cy. Here a are now 3D unit vectors. The 
mean field Hamiltonian is again 



TO 



- Jto ■ a — eff ■ a, 



(53) 
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where m is the magnetization, and this time we 
parametrize a as (sin 9 cos (j), cos 6 sin cos 0) . Therefore, 
the mean field equation reads 



= AVr, 



J a exp((3Jrn ■ a + Perjay) sin 9d0d(j) 
J exp(/3 Jm • (7 + Perjay) sin 0d0d<p 

where we have chosen 

m = (to sin ip cos x, w sin ip sin cos ip) . 



FH{m), 
(54) 



Consider first the case, when e = 0. By symmetry, the 
solutions of the mean field equation in this case form a 
sphere for /3 > 

Suppose that the radius of the sphere is rriQ . To find 
f3^^ analyticaUy, we can use an argument that is similar 
to the one that we have used in the case of the XY model. 
The existence of critical temperature gives the condition 



(55) 



which implies 



Pc J- 



The behavior of the magnetization as f3 approaches , 
is given by 



.o« = ^i(/3rr^(/3-/30^)i 



(56) 



Note that in the case of the XY model, we also found 
a similar behavior of the magnetization near its critical 
temperature. 

In the presence of disorder, i.e. e ^ 0, we obtain the 
correction to the critical temperature as 



3 9 6 / 9\ 



(3: 



■OH 



5 J3 



(57) 



which is again qualitatively siumilar to the situation in 
the XY model. The magnetization near the critical tem- 
perature is also decreased by an order of e^. The relation 
between the magnetization without disorder (to^), and 
that in the presence of disorder (to^) is 



1 - 



1 



12J2(7 



(58) 



Acknowledgments 

We acknowledge support from the DFG (SFB 407, 
SPP1078 and SPP1116, 436POL), Spanish Ministerio de 
Ciencia y Tecnologia grants FlS-2005-04627 and "Ramon 
y Cajal" , Acciones Integradas and Consolider QOIT, the 
ESP Program QUDEDIS, and EU IP SCALA. J. W. was 
partially supported by the NSF grant DMS 0623941. 



^ P.W. Anderson, Basic Notions of Condensed Matter 
Physics (Westview Press, Colorado, 1984); P. A. Lee and 
R.V. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985); R. 
Zallen, The physics of amorphous solids (Wiley, New York, 
1998). 

^ For a discussion of the recent development in the area 
of disordered quantum gases, see e.g. V. Ahufinger, 
L.Sanchez-Palencia, A. Kantian, A. Sanpera, and M. 
Lewenstein, Phys. Rev. A 72, 063616 (2005); M. Lewen- 
stein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen(De), 
and U. Sen, Adv. in Phys. 56, 243 (2007). 

^ M. Mezard, G. Parisi, and M.A. Virasoro, Spin Glass The- 
ory and Beyond (World Scientific, Singapore, 1987); S. 
Sachdev, Quantum Phase Transitions (Cambridge Univer- 
sity Press, Cambridge, 1999); D. Chowdhury, Spin Glasses 
and other Frustrated Systems (Wiley, New York, 1986). 

* D.J. Amit, Modeling Brain Function (Cambridge Univer- 
sity Press, Cambridge, 1989). 

^ A. Aharony and D. Stauffer, Introduction to percolation 
Theory (Taylor & Francis, London, 1994); G. Grimmett, 
Percolation (Springer, Berlin, 1999). 

^ A. Auerbach, Interacting electrons and Quantum mag- 
netism (Springer, New York, 1994). 

P.W. Anderson, Phys. Rev. 109, 1492 (1958); E. Abra- 
hams, P.W. Anderson, D.C. Licciardello, and T.V. Ra- 
makrishnan, Phys. Rev. Lett. 42, 673 (1979). 



* Y. Nagaoka and H. Fukuyama (Eds.), Anderson Localiza- 
tion, Springer Series in Solid State Sciences 39, (Springer, 
Heidelberg, 1982); T. Ando and H. Fukuyama (Eds.), An- 
derson Localization, Springer Proceedings of Physics 28, 
(Springer, Heidelberg, 1988). 

^ Y. Imry and S. Ma, Phys. Rev. Lett. 35, 1399 (1975). 
° J.Z. Imbrie, Phys. Rev. Lett. 53, 1747 (1984); J. Bricmont 

and A. Kupiainen, ibid. 59, 1929 (1987). 
^ M. Aizenman and J. Wehr, Phys. Rev. Lett. 62, 2503 

(1989); M. Aizenman and J. Wehr, Comm. Math. Phys. 

130, 489 (1990). 
^ J. Wehr, A. Niederberger, L. Sanchez-Palencia, and M. 

Lewenstein, Phys. Rev. B 74, 224448 (2006). 
^ C.J. Thompson, Classical equilibrium statistical mechanics 

(Clarendon Press, Oxford, 1988). 

* A. Aharony, Phys. Rev. B 18, 3328 (1978); D.E. Feldman, 
J. Phys. A 31, L177 (1998); B.J. Minchau and R.A. Feb 
covits, Phys. Rev. B 32, 3081 (1985). 

^ D.A. Abanin, P.A. Lee, and L.S. Levitov Phys. Rev. Lett. 
98, 156801 (2007). 

LA. Fomin, J. Low Temp. Phys. 134, 97 (2005); JETP 
Lett. 85, 434 (2007); see also C.E. Volovik, JETP Lett. 
81, 647 (2005). 

^ A. Niederberger, T. Schulte, J. Wehr, M. Lewenstein, L. 

Sanchez-Palencia, and K. Sacha, arXiv:07 07.0675 
® D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 



12 



(1966); P.C. Hohenbcrg, Phys. Rev. 158, 383 (1967). 
J. Zinri- Justin, Quantum Field Theory and Critical Phe- 
nomena, (Oxford Science Publication, Oxford, 1989). 
J. Frohlich, B. Simon, and T. Spencer, Comm. Math. Phys. 
50, 79 (1976); For more general proofs, see T. Balaban, 



Comm. Math. Phys. 167, 103 (1995); Comm. Math. Phys. 
182, 675 (1996). 

Handbook of Mathematical Functions, cds. M. Abramowitz 
and LA. Stegun (Dover, New York, 1970). 



